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We present exact and explicit form of the kernels K(x, y) appearing in the theory of energy 
correlations in the ensembles of Hermitian random matrices with Gaussian probability distribution, 
see E. Brezin and S. Hikami, Phys. Rev. E 57, 4140 and E 58, 7176 (1998). In obtaining this result 
we have exploited the analogy with the method of producing exact forms of two-sided, symmetric 
Levy stable laws, presented by us recently. This result is valid for arbitrary values of parameters in 
question. We furnish analytical and graphical representations of physical quantities calculated from 
K(x,y)'s. 

PACS numbers: 05.40.-a, 05.40.Fb 
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I. INTRODUCTION 

We consider in this work certain extensions of results 
obtained by Brezin and Hikami in hereafter referred 
to as BH. These authors have developed a theory of en- 
ergy correlations in the ensembles of Hermitian random 
matrices characterized by Gaussian probability distribu- 
tion. In the following we adopt the notation used by BH. 
In Sec. Ill of BH the Hamiltonian has a deterministic 
part Hq with the spectrum consisting of two distinct de- 
generated levels and the random part V. This produces 
the singularity in the density of states p(x) at the origin. 
The basic results of Sec. Ill of BH are contained in the 
form of the kernel K(x,y), see Eq. (3.21) of BH: 



K(x,y) 
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where the functions 
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(x) and ip{x) are defined by 
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where in Eq. ([2]) 



ei, see Eqs. (3.7) and (3.13) 



of BH. Evidently, the functions <p(x) and ij)(x) satisfy 
4>{—x) = cf>(x) and ^(—x) = —ijj(x). The BH provide, for 
both (j)(x) and ip(x), Taylor expansions for small x and 
the asymptotic behaviour at large x, see [l|, 0. These 



approximations are subsequently used to calculate the 
density of states given by p(x) = K(x,x) as 



p(x) 



4>'{x)^"{x) - 4>"(x)$'(x)+x4>(x)ijj(x) 



(3) 

as well as the connected correlation function for two 
eigenvalues symmetric with respect to the origin, given 
by 



p c (x,-x) = p c (x) 



K(x, —x) 



(4) 



Here we furnish exact formulas for <p(x) and ip( x ) thus 
enabling to visualize graphically Eqs. <j3j> and (|4|) for all 
the values of parameters. The exact forms of <p(x) and 
ij)(x) are 
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In Eqs. §5§ and © 0F2 is the generalized hypergeometric 

function of type p F q (jjj P ^ z^j as defined, for example in 

with (a p ) = ai, «2, ■ • ■ , ct p etc. 
In Sec. V of BH more involved conditions on Hq are 
chosen and this yields modified form of functions <j)(x) 
and ij){x) given by Eqs. (5.7) and (5.8) of BH respectively 
as 
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where in Eq. ((5J the integration path C consists of four 
lines of steepest descent starting from the origin to in- 

finity. On these lines t is replaced by t = e 3 1' and 

t = e 1 3 t', respectively. We have evaluated the inte- 
grals in Eqs. ([7]) and (JSJ exactly and their explicit forms 
will be reported in the part B of our Sec. II below. 

The paper is structured as follows: in Sec. II we derive 
the main results for the kernel K(x, y) and for its rela- 
tive K(x,y) built on 4>(x) and ip(x), an d summarize the 
asymptotics of <f){x), %j){x), <j)(x) and ip(x). In Sec. Ill 
we give graphical representations of these functions, of 
the both kernels and of their associated quantities /5(x), 
p(x), p c (x) and p c {x). We compare also graphically the 
exact results and the asymptotics of Sec. II. The Sec. IV 
contains discussion and conclusions. 
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where < Re(s) < 1. In Eq. (jllj) we have used the 
definition of gamma function along with the formulas 
2.3.3.1 and 2.3.18.2 in [ll|. We employ now in Re[p£(s)] 
the second Euler's reflection formula T\h + z)F(^ — 



, which leads, via Eq. (fTU|) , to 



Re[^(s)] = 4>{x) 



II. DERIVATION OF THE MAIN RESULTS 

A. We proceed to a streamlined derivation of the for- 
mulas (0 and ©• Thus we refer here to Sec. Ill of BH. 
The method is borrowed from our recent exposition of ex- 
act and explicit forms of one-sided [1] and two-sided [!, 0] 
Levy stable probability distribution functions. It is based 
on a combined use of Mcllin and Laplace transforms (for 
one-sided case) and of Mellin and Fourier transforms (for 
two-sided case). It translates directly to the case when 
the sought for functions are still of Levy-stable character 
but are not anymore probability distribution functions. 
According to [7| here the condition of positivity is not ful- 
filled and 4>(x) and ip( x ) are signed oscillatory functions. 
Utility of these functions in a number of applications has 
been strongly emphasized in Q . 

We define a complex function 



94(x) = - 



' 4 e" 



'dy. 
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It is clear that Rc[g4(x)] = 4>{x) and Im[g4(x)] is related 
to V»(x). We define the Mellin transform of a function 
f(x), for s complex, as 



M[f(x);s}=r(s) = / x s - 1 f(x)ds, 
Jo 

and its inverse as 



M-^ris^x] = f{x) = — / x- s r(s)ds, (10) 
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see [Io| for the appropriate choice of integration contour 
L in Eq. (fTTJ|) . 

In what follows we calculate the Mellin transform of 
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After changing variable u = Eq. (fl~3|) becomes 
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where in obtaining Eq. (|14[) we have used in the integrand 
of Eq. (|13|) three times the Gauss-Lcgcndrc multiplication 
formula for gamma functions. The form of Eq. (|14|) per- 
mits to naturally encode it as a certain Mcijcr's G func- 
tionjp, for which we use the conversion formula 16.17.2 
of jl3| in terms of generalized hypergeometric functions 
p F q (see also [l2[) and the final expression in the form of 
Eq. results. In order to meaningfully define the Mei- 
jcr's G function from Eqs. (|13j) and (fl"4")) the integration 
contours must be specified. This applies also to the con- 
tours in Eqs. (fT5 ]l . (flB ]) . fi ty and ® below. See HG3 
for a clear exposition of this procedure. 
The calculation of ip(x) follows a similar pattern. 

In Im[<74(s)] we apply the first Euler's reflexion formula 
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inversion of Eq. (fTU|) 



and it yields, after the Mcllin 
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Eq. ((T5|) will be further transformed with the change of 
variable u = and subsequent three-fold application of 
Gauss-Legendrc multiplication formula in the integrand 
of Eq. (|15|) , to arrive at 



Im[ 54 (y)] = 
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which is again clearly evocative of appropriate Meijer's G 
function. We apply anew the conversion formula 16.17.2 
of [13] to it and this yields: 



Im[.94(y)] 
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The function ^(x) will now be related to Eq. (JTTJ) . To this 
end we rewrite Eq. fl2J) using sinh(z) = — e 2 sin (e l 2 z) 
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which gives, according to Eq. (jTTj) 
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B. We refer now to the Sec. V of BH and consider in 
more detail our Eqs. ([7]) and (J5J. 

In close analogy to Eq. ([S]) we define 

1 r°° y 6 
9e(x) = — / e-Te-^dy, (18) 



1 3Te- i Tr( s )r(if) 1 



whose Mcllin transform turns out to be equal 

67T 

where < Rc(s) < 1. We repeat the reasoning applied 
to the gl(s) in Eq. (fTTjl and obtain (compare Eqs. (fTTj) 
and (|12| above): 
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Likewise, 
Im[g 6 (x)] 
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In order to represent Eq. (jHJ) as an integral over posi- 
tive axis we sum up four contributions as specified after 
Eq. (J5J and end up with 
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that gives with the help of Eq. (|20]l the final result 
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Thus, Eq. (5.6) of BH which is (we have removed the 
"hat" to be consistent with Eqs. © and ©) 



K(x,y) = 



^\x)^{y)~^\xW{y) + r{x) 

(22) 



x-y 

x 0'(x)v (3) (y) + 0(^^ (4) (y) 

gives the density of states 

p(x) = -x</»(x)V(x)-2U^(x)V'(x)-0'(x)V' (4) (x) 



0"(x)V (3) (x)-0 (3 )(x)V'"(x) 



(23) 



where the first term is the previous expression stems 
from the exact differential equation for </>(x): cf>^(x) = 
-§</>(x). With </>(x) of Eq. {SJ and ^(x) of Eq. (ETJ the 
density of states p(x) and the quantity 

Pc (x) = -[A'(x,-x)] 2 (24) 

analogous to Eq. Q can be readily calculated. 
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C. Although exact forms of 4>(x), ijj(x), <j>{x) and ip( x ) 
render possible all sorts of calculations, it is still instruc- 
tive to consider the asymptotics of these functions. We 
shall briefly comment on it now. In doing so it is practi- 
cal to introduce the so-called Levy functions of index a, 
a > 0, which are two-sided symmetric Levy stable dis- 
tributions denoted in [H, @ as g(a,0;x), and considered 
there under the condition < a < 2. They are spe- 
cial case of general two-sided Levy function g(a,/3;x), 
f3 — a, a > 0. In the present work /? is set to zero and 
therefore 0, i, 



g(a,0;x) 



e 1 cos(xt)dt, 



where a > and — oo < x < oo. As mentioned ear- 
lier only for < a < 2 are these functions probability 
distribution functions (pdf). For a > 2, g(a,0;x) are 
not pdf's and are signed oscillatory functions. One sees 
rapidly that 



and 



fa)=2 1 ' 3 g(4,0;2 1/2 x), 



<f>(x) = 3 1/6 .g(6,0;3 1/6 x). 



Although initially dismissed as not being the pdf's, the 
functions g(a,0;x) reappear in the literature again. A 
considerable amount of information is contained in [H, . 
In 0, we have obtained exact and explicit forms of 
g(a,0;x) for admissible rational a. which include a — 4 
and 6. We emphasize that our formulas for g(a, 0; x) fur- 
nish exact expressions even when they are not pdf's. We 
employ now Eq. (5.18) of to confirm the asymptotic 
behaviour for x — > oo in Eq. (3.10) of BH and (2.15) of 
0, that is 
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Eq. (5.21) of gives the asymptotics of 4>(x) in the form 



<P(x) ~ (2/ ^l 5 \ l cxp[- K {x)] + exp [cos (%)k(x)] 
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x cos 



(26) 



for x -> oo with k(x) = (5/3) (x/2) 6/5 . 

Using the Taylor expansion of cosine in Re^^)] in 
Eq. (|18[) leads to the asymptotic behaviour of <j){x) at 
x — > 0, namely 



5 oo 1 

3~6 (-33) 



2tt 



r=0 



2/- 



(27) 



To find the asymptotic behaviour of tp(x) f° r x — > oo 
we have a ready-to- use formula 6.46 on p. 96 of [l4j], i.e. 
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which is valid for n = 2, 3, 4, . . ., arg(z) = or for n = 1, 

|arg(z)| < 2n, with A = Efci c fc, £ 
That gives 



i/2-X 
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Vi(x) 



yjz-i exp (fx 4 / 3 ) cos(^ 4 / 3 + f ), (28) 



which apart the factor 2, agrees with Eq. (3.16) of BH, 
and Eq. (2.21) of 0. 

By similar means we conjecture the asymptotics of ij}{x), 
x — > oo, to be 
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x sin 



|^ exp[-/c(ar)] - exp [-cos (^)k(x)] 
sin(f )k(x) + E]} + ... . (29) 



The asymptotic behaviour of <fi(x) at x — > can be ob- 
tained by taking Taylor expansion of Eq. (|2ip . That 
yields to 
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Substituting the asymptotics at x — > oo of 0(x), 4>{x) 
and 0(cc), V^) it is easy to estimate the behaviour of p(x) 
(see Eq. ©) and (see Eq. ([23])) for large cc. That 
gives the expressions whose first terms with the leading 
powers are proportional to x 1 / 3 and x 1 / 5 , for p(x) and 
p(x), respectively. 

In the following Section we shall compare this asymp- 
totics with the exact results. 



III. GRAPHICAL REPRESENTATION 

As the formulas from Sees. II, A and II, B arc valid for 
arbitrary x, we have carried out graphical presentation 
of selected quantities. In Fig. [T]we display the oscillating 
functions <j>{x) and 4>{x), see Eq. ([5]) and (TT9")) respectively. 
The amplitude of these oscillators dies down very rapidly 
with x. In Fig. [2] we present the oscillating functions 
-0(x) and ip(x), see Eq. and (f21~j) . respectively. Here, 
the amplitudes increase drastically as a functions of x. 
In Fig. |3] some sections through the surface K(x, y) (see 
Eq. (JTJ) ) arc presented, for fixed values of y, as a function 
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of x. We have also obtained similar curves for K(x,y) 
from Eq. (|2"2"j) but we shall not present them here. In 
Fig. |4]we have presented the densities of states p(x) and 
p(x) from Eq. ([3]) and (j2"5|) respectively. They present 
a non-monotonic behaviour with characteristic pockets 
where local minima appear. In Fig. [5] the same quanti- 
ties are compared with their asymptotics for large x, see 
Sec. II, C. For transparency, in Fig. [5] we have retained 
for the dotted curves only the leading power-law terms. 
The oscillatory contributions to the asymptotics are not 
displayed. We stress however that the oscillations in the 
exact solutions (see continuous curves) agree very well 
with the oscillations given by asymptotics, especially for 
x > 20. In Fig. [6] the correlation function —p c (x) is dis- 
played, see Eq. (j4j. In Fig. [7] we illustrate the correlation 
function —p c (x), see Eq. (f!H)) . 




FIG. 1: (Color online) Plot of the functions 4>(x), see Eq. |0, 
red line, and <f>(x), see Eq. (JT5J, blue line. 
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FIG. 2: (Color online) Plot of the functions 4>(x), see Eq. ((6J, 
red line, and ip(x), see Eq. (|21[) . blue line. 
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IV. DISCUSSION AND CONCLUSIONS 

The key results of present investigation arc contained 
in the exact expressions for the kernel K{x, y) in Eq. ([1]) 
and for the kernel K(x,y) in Eq. (f22j) . This allows one to 
obtain with an ease the density of states p{x) and p(x) 
which display a characteristic oscillatory non-monotonic 
behaviour, with the infinite scries of local minima. 

One observes that both pairs of relevant functions 
0(2;), ip{x) and </>(x), ip{x) can be expressed with a par- 
ticular subset of generalized hypcrgcomctric functions, 

namely the type o^n ( ( c _ ) ' n ~ 1 3 2, . . ., where c„ are 
positive parameters. These functions have been investi- 
gated in the past [l4], [13 an d are sometimes referred to 



FIG. 3: (Color online) Some sections through the surface 
K(x,y), see Eq. |T]), as a function of x, for: y — 1/2 (red 
line), y = 1 (blue line), y = 3/2 (green line), y = 2 (violet 
line), y = 5/2 (magenta line). 



as hyper-Bessel functions |14| . They play a prominent 
role in the fractional calculus fl6j which in turn is a ba- 
sic tool to investigate non-standard diffusion in statistical 
mechanics (l7j . 

The functions encountered in the course of our evalu- 
ations can be shown to be related to oscillatory integrals 
used to describe wave propagation and optical diffrac- 
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FIG. 4: (Color online) Plot of densities of states: p(x), see 
Eq. Q, red line, and p(x), see Eq. ((231) . blue line. 
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FIG. 5: (Color online) Plot of densities of states: p(x) see 
Eq. 0, blue line, and its large- a; asymptotics, ~ 0.276a; 1 / 3 , 
dotted red line; p(x) see Eq. (|23[). black line and its large- 
x asymptotics, ~ 0.270a; 1 '' 5 , dotted brown line. For trans- 
parency the dotted curves represent the leading power ap- 
proximations, without the inclusion of the oscillating terms. 



tion. In particular, our function gAx) is closely related 
to the so-called Pearcey's integral [18[ . Its various forms 



and asymptotics can be found in Chap. 36 of |14[. 

We believe that the functions of Levy type will be also 
appearing in further studies of spectra of random matri- 
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FIG. 6: Plot of correlation function —p c (x), see Eq. (j4|). 
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FIG. 7: Plot of correlation function —p c (x), see Eq. ([24]) . 
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